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Abstract. We analyse a Far InfraRed (FIR) catalogue 
of galaxies at 60 /im with a flux limit of w 110 mJy, ex- 
tracted from a deep subsample of the IRAS Faint Source 
Survey. Monte-Carlo simulations and optical identification 
statistics are used to put constraints on the evolution of 
galaxies in the FIR. We find evidence for strong evolution 
of IRAS galaxies, in luminosity oc (1 4- z)^-^ ±o.2±o.3 q^, 
density oc (l-f 2)60 ±o.5±o.7 ^j. j^(60^m) > 150 mJy, in 
agreement with previous claims. An excess of rather red 
faint optical counterparts with B > 18, consistent with 
the above evolution, is detected. We interpret these re- 
sults as strong evolution at recent times among the star- 
burst (or dusty AGN) population of merging/interacting 
galaxies. Most of these objects at moderate redshifts may 
pass unnoticed among the population of massive spirals 
in broad-band optical surveys, because of large amounts 
of dust extinction in their central regions. A possible link 
with the strong evolution observed in the optical for blue 
sub-M* galaxies is discussed. 

Key words: Extragalactic astronomy: infrared: galaxies 
- galaxies: evolution - galaxies: starburst - galaxies: in- 
teractions ~ galaxies: fundamental parameters 



1. Introduction 

This is the second paper in a series where we try pho- 
tometricaly to quantify galaxy evolution at low redshifts 
(z < 0.3), from reasonably large and homogeneous sam- 
ples. We have shown in Bertin & Dennefeld ( |1996 , paper 
I) that optical galaxy counts were consistent with little 
evolution down to Bj = 21, corresponding to 2: sa 0.2 for 
M* galaxies. In this paper, we investigate the case of the 
far-infrared domain. 



The Far-InfraRed (FIR) is very sensitive to ongoing 
star-formation activity in galaxies. Warm dust, essentially 
heated by the blue/UV continuum of young OB stars, 
makes a large contribution to the total FIR output of late- 
type objects. At the same time, a significant heating by 
an Active Galactic Nucleus (AGN) is possible in the most 
luminous FIR-emitting objects, although the respective 
contributions of star-formation and AGNs are not yet well 
determined. In this context, a statistical analysis of a deep 
and homogeneous FIR-selected galaxy sample should eas- 
ily allow one to trace evolution of star-formation and/or 
AGN activity with redshift. 

The IRAS (InfraRed Astronomical Satellite) survey, 
with its four infrared channels at 12, 25, 60 and 100 iim 
provides a unique tool for such studies, despite its moder- 
ate sensitivity. Based on the Point Source Catalog (PSC) 
and the more recent Faint Source Catalog (FSC) extracted 
from the IRAS data, several workers have investigated 
the evolution hypothesis on subsamples reaching various 
depths. The situation is still somewhat unclear, although 
evolution seems to be detected at 60 fj,m in most cases. 
Using PSC data, with solid identifications and radial ve- 



locities, Saunders et al. (1990) claimed a very strong evo- 
lution effect in their sample, while Fisher et al. (1992) 
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found theirs compatible with no evolution (the latter has 
however a somewhat lower media n reds hift). Using the 
deeper FSC data, Lonsdale et al. ( 1990 ) found evidence 
for strong evolution from source number counts. Oliver 
et al. ( 1995 ) reported also evolution in their spectroscopic 
subsample, though at a slightly milder rate. Two deeper 
studies (Hacking & Houck 1987 and Gregorich et al. 1995 ) 
made on much smaller specific areas (a few tens of square 
degrees) suggest a large excess of faint detections. 

In order to clarify the problem and to complete these 
studies by providing the "missing link" between exten- 
sive, shaUow catalogs (/^(60/^m) > 0.2 Jy) and small, very 
deep studies (/,y(60/im) ^ 100 mJy), we have undertaken 
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an identification program of IRAS sources at 60 /im on 
some selected areas in tfie sky, focusing ourselves espe- 
cially on the faintest flux domain achievable with the Faint 
Source Survey (FSS) "plates": /^(60Aim) « 120 mJy with 
a signal-to-noise ratio of 4. Our purpose is twofold: (1) to 
estimate the reliability of the IRAS survey at its very faint 
end, and (2), to characterize the population of IRAS galax- 
ies in the 100 mJy domain (where evolution is expected 
to begin to show up in models of FIR galaxy counts). 

At this low S /N level, the data cannot be interpreted in 
a straightforward way: completeness, reliability and mea- 
surement biases must be quantified accurately. We rely 
here on both Monte-Carlo simulations and optical identi- 
fications to estimate these effects. 

The paper is organized as follows. In s|| we describe 
the selection of the infrared "real" sample, sjs] is devoted 
to the description of our Monte-Carlo simulations of IRAS 
catalogs. These are then used in §^ to interpret statistics 
on the infrared data (confusion noise, number counts). In 
^ we rely on optical identifications to put further con- 
straints on the nature of faint sources. We summarize our 
results in §^ and compare them to those of other studies. 
In ^ we finally discuss the implications in the framework 
of galaxy evolution. 

2. The data 

The data in this study come from the Very Faint Source 
Sample (VFSS). The VFSS was generated by one of us 
(MM) after the completion of the FSS. The data in the 
FSS were searched for regions where 1) The IR cirrus indi- 
cators NOISRAT and NOISCOR (Moshir et al. 1992) were 
small, 2) the spacecraft coverage was at least 4 HCONs 
and 3) there was a minimum number of nearby galax- 
ies. These conditions led to a list of FSS plates which 
were then processed through the normal processing pipe- 
line. Their average "instrumental" noise is 25 mJy instead 
of about 40 mJy for more typical FSS plates. Since with 
these conditions, the reliability of the sample had been 
verified to be high through spot checks, the extraction 
thresholds for the normal FSS pipeline (which were set 
conservatively) were reduced (from S/N of 3 to 2.5) in or- 
der to increase the completeness of the VFSS down to 110 
mJy at 60 /xm (and similar increases in completeness in 
the other 3 bands). The resulting VFSS database covers 
400 square degrees in 4 separate contiguous regions (Table 
|l|) and contains approximately 2,500 sources detected at 
60 /xm. The VFSS is about 2 times deeper than the FSC 
and contains many unique sources since they can not be 
found in any other IRAS catalogue or databaseEl. 

Although we have only upper limits at 12, 25 and 
100 /im for the great majority of the 60 /im detections. 



Table 1. IRAS fields 



Field a (2000.0) 5 (2000.0) Area 





h 


m 




sq. degrees) 


Nl 


09 


20 


+78 20 


80.8 


N2 


15 


25 


+55 00 


113.5 


N3 


17 


50 


+47 30 


83.2 


S 


03 


50 


-47 00 


133.7 



a few stars lie in the sample. We removed them by re- 
questing infrared sources to have /i/(25/im) < /,y(60/im). 
This simple criterion proves to filter efficiently most galac- 
tic objects, but might also discard a few Seyfert galaxies 
(e.g. Dennefeld & Veron 1986). However, as only 2% of 
our detections do not pass the filter criterion (all identi- 
fied as bright stellar sources in our optical subsample), we 
consider this effect to be negligible. 



^ Although the field studied by Hacking & Houck (1987) has 
one of the highest spacecraft coverages and thus passes condi- 
tion 2, it fails to pass condition 1 for cirrus contamination and 
was automatically excluded. 



3. The "Monte-Carlo" 60 fim sample 

About 2/3 of our detections at 60 /xm have a S/N lower 
than 5, i.e. a photometric accuracy worse than 20%. At 
these levels, one expects the data to be seriously affected 
by completeness and reliability problems. Important pho- 
tometric errors (if they are not purely fractional) do sig- 
nificantly bias number counts (Eddington 1913), as well 
as many other statistics that can be drawn from a flux 
limited survey. Murdoch et al. ( 1973[ ) have tabulated the 
Eddington bias, in the typical case of an underlying Eu- 
clidean slope and Gaussian flux errors, down to a "safe" 
5cr level. For detections lying below this lev el, the sim- 
ple analytical correction used by Oliver et al. ( 1995 ) from 
the Murdoch et al. calculations, cannot be applied. The 
only escape is to proceed through Monte-Carlo simula- 
tions. This can be easily done here as both the instru- 
mental and total noise levels are nearly constant over the 
surveyed zone. We have therefore created mock infrared 
images mimicking FSS plates, on which was run a source 
extraction algorithm similar to the original one used for 
the FSC. Statistics drawn from the simulation catalogs 
could then be directly compared with those from the real 
observations, with the hope of discriminating against dif- 
ferent models of the IRAS galaxy population. 

The Monte-Carlo approach has some limitations at 
these low S/N levels, mainly because assumptions have 
to be made about the noise distribution. For this rea- 
son, we included only two components in the simulations: 
"instrumental noise" , which we considered to be Gaus- 
sian, and galaxies. Given the very high coverage in our 
fields, the data conditioning involved in the production of 
FSS plates (including de-glitching, median filtering and 
trimmed average) should have removed all significantly 
non-Gaussian features from the instrumental noise distri- 
bution. Filaments from infrared cirrus (Low et al. 1984) 



are more of a threat, and can easily create false detections 
by themselves. In what follows, we will keep in mind the 
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possible existence of such noise components. But one can 
remark that the effect of transient noise sources like cir- 
rus peaks and glitches, whose distributions are strongly 
skewed toward positive values, is essentially to lower the 
reliability of detections. 

3.1. The galaxy population 

We modeled the flux distribution of infrared galaxies as- 
suming a c onstan t comoving number density, and a Saun- 
ders et al. (199C) Luminosity Function (LF): 



0(i6o) = C 



exp 



1 

'2^ 



logio 



L 

L* 



(1) 



where Lgo = i^L^ at 60 /xm, a = 1.09, a = 0.724 
and L* — 2.9510® H'^Lq. These best-fit parameters de- 
pend slightly on the value of Hq, which was conveniently 
set to GOkm.s^^.Mpc""'^ in our simulation, close to the 
value of 66km.s~^.Mpc~^ adopted by Saunders et al.. 
(j){Leo) was normalized so as to fit the QDOT galaxy 
counts (Rowan-Robinson et al. 1991) over the domain 
1 < /i/(60/im) < 2 Jy. It therefore assumes that evolu- 
tionary effects are negligible at this flux level. We adopted 



the luminosity-dependent fit of Hacking & Houck (1987) 
to the Soifer et al. (1987) Bright Galaxy Sample, for the 
distribution of spectral indices. Pure luminosity evolution 
was introduced in some of the realizations by making the 
L* parameter evolve as (1 -|- z)'^ . This simple form is not 
introduced here with strong physical justification, but only 
to provide a convenient estimator for the evolutionary rate 
of the galaxy population in the infrared. Luminosity evo- 
lution was preferred to density evolution as it is more con- 
venient to simulate through Monte-Carlo methods. At the 
depth reached by our data, significant evolution can only 
be probed for FIR-luminous galaxies. As for these galaxies, 
the LF behaves much like a power-law (Eq. |l|), very similar 
results are obtained for pure density evolution. We gave 
to the "cool" and "warm" components (Saunders et al. 



1990) the same rate of evolution for simplicity. In fact the 
"warm" (starburst) component proves to be the essential 
contributor to the part of the LF where evolution can be 
probed here. 

The spatial distribution of galaxies was simulated us- 
ing a three-dimensional log-normal density field (Coles & 
Jones 1991), and the Saunders et al. (1992) autocorrela- 
tion function truncated at 20 Mpc and kept constant in 
spatial coordinates. The log-normal model gives a honest 
fit to the count probabilities in the 1.2 Jy sample (Bouchet 
et al. 1993) and, although it may be unable to reproduce 



very high density regions of features like sheets and fila- 
ments, it provides an interesting alternative to other sim- 
ple "to y univ erse" construction methods (e.g. Soneira & 
Peebles [l978l Benn & Wall 1995). 

The overall simulated galaxy sample was then 
placed in an Einstein-de Sitter universe, within a rod 
150 h^^ Mpc wide, extending up to z = 2. This redshift 



limit was chosen here for practical purpose. It does not 
affect by more than a few % counts with /i,(60/im) > 
10 mJy for the highest evolution rates considered in this 
paper. Although (projected) clustering was obvious on 
simulation images as seen by eye, none of the statisti- 
cal measurements discussed in this paper proved to be 
affected significantly if we switched to a simple Pois- 
sonian distribution. Nevertheless, all numbers presented 
here were derived from simulations with clustering. 

3.2. The images 

For the sake of simplicity, we simulated directly IRAS im- 
ages by projecting galaxies over the frame, without in- 
cluding the complex co-adding of individual scans. The 
simulated images were given the same pixel size as the 
original FSS plates: 0.5'. This generously samples the re- 
constructed IRAS "point spread function" , thereby re- 
moving the need for oversampling when projecting the 
infrared sources before convolution. At the same stage, 
we added white noise to account for instrumental noise. 
Images were convolved by Fast-Fourier Transform with a 
3.7' X 2.1 'Gaussian beam (a reasonable fit to the 60 fim 
FSS templates at these ecliptic latitudes), with the minor 
axis (average "in-scan" direction) aligned with the y-axis. 
Finally, we suppressed the mean background and low spa- 
tial frequencies by applying the standard 60 /im zero-sum 
median filter of the FSS (Moshir et al. |1992| ) along the 
y-axis of the image. In the FSS, this point-source filtering 
was normally done on individual IRAS scans, which im- 
plies some angular dispersion in the final images at high 
ecliptic latitude. Anyway, our simulated "in-scan" point- 
source profiles are in excellent agreement with those from 
the original plates. 

3.3. Source extraction 

The mock catalogue was created from the simulated im- 
ages by running a slightly modified version of the SExtrac- 
tor (Bertin & Arnouts 1996) photometry program, tuned 
to mimic closely the original FSC "Extractor" from IPAC 
(Moshir et al. 1992). In particular, the standard "clipped 
rms" estimate of background noise in SExtractor was re- 
placed by the "68%-quantile" from the FSC extractor, and 
the integrated flux of detected sources by the peak pixel 
intensity, as for the real data. 

4. Infrared statistics 

4-.1. Confusion noise 

It is interesting to compare the average 68%-quantile mea- 
sured on simulations with those from the real flux-maps 
(Table ^ . This parameter is particularly sensitive to the 
density of undetected sources: in quadrature, about 2/3 
of the confusion noise is expected to come from sources 
in the 10-100 mJy flux range with the models considered 
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here. The 68%-quantile is fairly stable in our simulations 
and in the real measurements, and enables us to constrain 
the number of these faint galaxies. Model-dependent con- 
clusions about the evolution rate can be drawn from it 
(there is of course a dependence on model normalisa- 
tion itself, but number counts of bright sources should 
restrict this uncertainty to less than « 10%). As Table |2] 
shows, the fluctuations predicted by a model evolving as 
L* oc (1 -|- z)^, with Q ^ 4 are too large (the uncertainty 
quoted here represents the rms of the 68%-quantile over 
the southern subset area, where it reaches its minimum 
value). The Q = 3 and Q = 3.5 models agree well with 
the data if there is no other contribution than pure instru- 
mental noise affects in our FSS fields. Milder-evolution and 
no-evolution models may require additional noise, cirrus 
for example, to fit the data. 



Table 2. 68% quantile statistics in mjy (including 25 mjy rms 
instrumental noise). 



VFSS data 


no evol. model 


= 3 


= 3.5 


= 4 


30.2 ± 1.2" 


28.1 


30.1 


30.7 


32.8 



standard deviation (see text). 



Some very faint cirrus can unambiguously be seen 
in the unfiltered 60 /im maps, which means they may 
somehow contr ibute to the background confusion noise. 
Gaultier et al. ( 1992 ) have estimated the confusion noise 
at 60 /im due to cirrus to be about 7 mJy rms, at high 
galactic latitude in the FSC. This might be taken as a 
typical value for our catalog, as our fields are located at 
only |6| « 35, but are amongst the "cleanest" of the in- 
frared sky concerning cirrus emission, and are comparable 
to higher galactic latitude zones. Because the infrared sur- 
face brightness distribution of cirrus has generally much 
larger wings than a Gaussian, one should expect their con- 
tribution to the total 68%-quantile to be less than 7 mJy 
in quadrature. However, given the large uncertainties in 
the real importance of cirrus emission in our fields, we do 
not discard the possibility of significant cirrus contamina- 
tion of both background and individual source measure- 
ments. 

In any case, this fluctuation analysis is consistent with 
previous estimates that used the integrated infrared back- 



ground (Oliver et al. 1995), although it is slightly more 
stringent for the Q = 4, Zf = 2 model which is only 
marginally compatible with our data. 

4-. 2. Number counts 

Deriving (raw) infrared number counts is straightfor- 
ward here, as we are dealing with a homogeneous set of 
data. Comparison with our Monte-Carlo samples is shown 
in Fig. 1^. The Eddington bias is clearly visible below 
/y(60/xm) « 200 mJy. The best fit value for Q proves 



to depend on the flux range considered. For f^{60iim) > 
150 mJy, i.e., for sources above with S/N higher than 5, we 
find Q — 3.2 ±0.2, which however raises to Q = 3.6 ±0.3 if 
we consider only the 0.15-0.2 Jy range. In the more uncer- 
tain 0.11-0.15 Jy domain, we find Q = 3.9 ±0.2. Account- 
ing for a ±10% possible error in the model normalisation 
yields a « ±0.3 additional uncertainty to the values of 
Q quoted above. These numbers may indicate that the 
(1 + z)"^ dependence in luminosity is inappropriate for de- 
scribing the effect seen here. However, we shall remember 
that, as stressed earlier, assumptions about the noise dis- 
tribution become very critical in the 0.11-0.15 Jy bin and 
will affect conclusions made with low S/N data. 

This brings the question of whether our fields represent 
a "fair sample" of the Universe, i.e. is there a relative ex- 
cess of faint or bright sources due to large scale-structures 
which could distort the statistics (see e.g. Lonsdale & 



Hacking |1990| )? Fortunately this point can be addressed 
using our bright optical galaxy counts in the same areas 
(paper I). Expressed with the same magnitude scale, they 
are similar to those done on much larger catalogs (APM, 
COSMOS) in the 17 < Bj < 21 range, where most of the 
optical counterparts to faint IRAS sources lie. A strong 
distortion of galaxy counts induced by the presence of 
large scale structures is thus excluded. 

An interesting question is: what happens to number 
counts if we add some Gaussian noise to the no-evolution 
model, in order to reach the observed intensity of fluc- 
tuations? The answer is: almost nothing, as can be seen 
in Fig. 0. As expected, this test indicates that the no- 
evolution model could be saved only if the distribution 
of "hidden" noise components is extremely non-Gaussian, 
creating noise peaks without increasing too much the 68% 
quantile. 

5. Optical identifications 

In order to find out whether the faint source excess mea- 
sured in the previous sections is real or contaminated by 
noise peaks, some identification work is necessary. Sys- 
tematic studies conducted at higher IRAS flux limits (e.g. 
Wolstencroft et al. |l986l |l992l Sutherland & Saunders 



1992D have shown that the optical domain is well adapted 
to this task. 



5.1. The data 

Thirty percent (6 Schmidt plates) of the VFSS were anal- 
ysed using photographic data in the blue Bj and red Rp 
passbands. The sensitivity of deep survey plates is well 
matched to our IRAS study, as the VFSS error box at the 
faintest flux levels occupies roughly 1 square-arcminute on 
the sky, and contains already an average of « 0.3 galax- 
ies brighter than Bj = 21. Full details about the optical 
catalog can be found in Paper I; we summarize here its 
main features. This catalogue is quite small compared to 
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The PSC/QD OT (R owan- Robinson et al. 
(Lonsdale et al. 199C ) counts have been added for comparison. 
The FSC counts are not corrected for the Eddington bias, and 
so their two faintest bins are probably overestimated. For clar- 
ity, the Monte-Carlo models were plotted only for fluxes below 
1 Jy. solid: no-evolution -I- Gaussian noise (see text); dashes: 
evolution with Q = 3; dots: Q — 3.5; long dashes: Q — 4. Error 
bars are Poissonian uncertainties. Monte-Carlo models become 
very noisy brightward of 1 Jy, and have not been drawn there 
for clarity. Note that the bin size is not constant. 



large galaxy surveys like APM (e.g. Maddox et al. 1990) 
or COSMOS (e.g. Heydon-Dumbleton et al. |1989[ ), but is 
well calibrated (within « 0.1 mag, thanks to a high density 
of photometric standards). The completeness for galaxies 
is expected to be better than a conservative 90% at least 
for 17 < Bj < 20. Brighter than this, there might be some 
loss due to inaccurate star-galaxy separation, and under- 
estimation of fluxes (w 0.1 mag) because of photographic 
saturation effects. Positions are accurate within a fraction 
of arc-second. This is sufficient for identifying faint IRAS 
sources, which have typical positional uncertainties of a 
fraction of arcminute. 

5.2. The simulated optical sample 

As for the infrared statistics presented in the previous sec- 
tions, extreme caution should be taken in interpreting the 
results of optical identifications. We don't know a priori 
what fraction of real IRAS sources can be identified, espe- 
cially given the large errors in /j,(60/ini) at our flux level. 
Here again, we addressed this question using a Monte- 
Carlo method, by simulating a sample of identified optical 
counterparts to the artificial sample described in 

5.2.1. Optical luminosity 

We first assigned some random optical luminosity to each 
infrared "Monte-Carlo" galaxy with theoretical flux re- 
ceived on Earth /i,(60/Lim) > 10 mJy, according to a suit- 
able optical/60 /im distribution. This is an easier task 



than it may appear, as IRAS and bright optical galaxy 
catalogs do have an overlap large enough so that one can 
build a local, bivariate luminosity function. The most com- 
prehensive study in that field was done by Saunders et al. 



(1990), who proposed a Gaussian distribution of the blue 
absolute luminosity, with standard deviation ctm = 0.72 
and mean 



HM = -19.9+5 log /i-f 0.3 ( log-^ - 11.24- log /i-^ ) .(2) 



Note that the original expression contains a sign error. 
M normally refers here to a Bj, magnitude, which cor- 
responds almost to our Bj magnitude system (Paper I) 
as we drop the average correction for internal extinction 
Ai = 0.21 (Rowan-Robinson [l987| ) . Although (|) is based 
on limited statistics and rather inaccurate optical mag- 
nitudes, we believe it does provide a globally correct de- 
scription of the link between the 60 /im and blue luminos- 
ity functions because: 1) relation (||) correctly reproduces 
both the L Pin/ Lb distribution for "normal" luminosity 
sources (Rowan- Robinson et al. |1987[ Bothun et al. p_989 ) 
and the "saturation" of optical luminosities observed in 
Ultraluminous InfraRed Galaxies (ULIRGs) (Smith et al. 
1987, Soifer et al. 1987); 2) integrating over all possible 
infrared luminosities leads to a blue luminosity function 
in excellent agreement with recent, local dctcrminationso 
(e.g. Ellis et al. 1996). The bivariate luminosity function 
based on (|]) and (||) is shown in Fig. ^ We assume rela- 
tion (H) to hold out to z « 0.3, corresponding to the most 
distant, luminous galaxies we might detect in numbers 
both in FIR and in the optical. This obviously supposes 
a negligible evolution of their dust content over the cor- 
responding period of time. We shall tackle this point in 



5.2.2. fc-corrections 

In order to give the right optical fc-corrections to our sim- 
ulated IRAS sources, one has to make some general as- 
sumption about their optical/UV spectral energy distri- 
bution (SED). The large majority of IRAS galaxies con- 
tained in FIR flux-limited samples are spirals (e.g. de 
Jong et al. |1984), w ith average morpho logica l type Sbc 



(Bothun et al. |l989i Sauvage & Thuan [19941 ). Although 
their rest-frame optical colours span a considerable range, 
it is difficult to flnd an obvious trend with FIR luminos- 



ity (Leech et al. |1988| , |l989i Hutchings & Neff |l99lD . An 
Sbc SED also proves to reproduce satisfactorily the ob- 
served optical continuum of 1.4-GHz sub-mJy starburst 
galaxies with B < 22 (Benn et al. 1993). Given the 
tight correlation existing between radio-continuum and 
far-infrared galaxy luminosities (de Jong et al. 1985, Helou 



^ It is interesting to note that this leads to a "high" normali- 
sation of the local blue luminosity function, in agreement with 
our bright optical counts (Paper I). 
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Fig. 2. Contour plot of the assumed bivariate 60 /im/Bj 
galaxy luminosit y fun ction, from the parameterization work of 
Saunders et al. ( |l990| ). There is a factor 2 in density from one 
contour to another. The luminosity range analysed by Saun- 
ders et al. is approximately delineated by the heavier contour. 
The dotted curve represents Eq. (^. Dashes indicate lines of 
constant Lqo/Lb- 



et al. |1985D , /^(1.4GHz) w 10~'^f^{60fim), these sources 
are representative of the faintest ones we might detect in 
our 60 /im sampl e. We have therefore adopted the Sbc 
SED from Pence (|l976|) to compute optical fc-corrections 
through our Bj and Kp passbands, and the average Sbc 
(Bj — Rp) = 0.96 typical rest-frame colour from Paper I. 
We will find in § 5^ other arguments justifying this choice. 

We finally end up with a large number of artificial 
sources having in addition to their "unbiased" /^(60/im) 
fluxes, Bj and Kp magnitudes that could directly go 
through the same complex IR detection and optical iden- 
tification processes, as real galaxies. 

5.3. Identification method 

Because of the low resolution of IRAS images 2' at 
60 /Lim), the reliability of the identification technique is 
crucial for finding the right optical counterparts to faint 
infrared sources. This is even more critical here, because 
of the particularly low infrared S/N and optical faint- 
ness for most of the candidates. We have adopted the 
likelihood-ratio method, inherited from radio-source iden- 
tifica tion p rograms (de Ruiter et al. 1977] , Prestage & Pea- 
cock 1983), which was successfully applied in previous 



IRAS studi es (W olstencroft et al. |l98e 
& Saunders |1992D. 



The positional uncertainty of an IRAS detection is es- 
sentially a function of the local reconstructed IRAS point 
spread function and the S/N. Assuming the error distribu- 
tion is roughly Gaussian, this defines the traditional error 
ellipses around the estimated object position, described by 
the ax, (Ty and 6\ respectively the rms error along the ma- 
jor axis (cross-scan), the rms error along the minor axis 
(in-scan), and position angle. Switching to the reduced 
coordinate: 



(3) 



we can express the probability for the true optical coun- 
terpart to lie between r and r -\- dr as 




Pr{r) dr ~ Q e^^ ' r dr 



(4) 



where Q is the probability for the infrared detection to 
exist and to be detectable in the optical catalogue. In or- 
der to check that this modeling fits the data correctly, 
we derived maximum-likelihood estimates of ax, ay and 
9 from optical identifications of the brightest infrared 
sources {fv{60 fj.m)> 0.2 Jy). We were able to recover the 
catalogue values within 10-20%. Of course, this does not 
mean that the modeling is adequate. There is indeed some 
evidence that the wings of the positional error distribution 
are larger than Gaussian ones (e.g. Sutherland & Saunders 
1992); we will quantify the consequences of this effect in 



the next section. As shown in Fig. ^, the Gaussian model 
proves to fit correctly the distribution of optically bright 
(Bj < 17) counterparts in our catalogue. One can how- 
ever already detect a slight excess (w 5% of the total) of 
identifications at about 3a along the cross-scan direction. 



1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


... 1 ... 1 


. . 1 . . . 




1992, Sutherland 





X/a, 



Fig. 3. Distribution of the 227 galaxies with Bj < 17 and 
r < 4 detected around infrared sources, along the major (left) 
and minor axes of the IRAS error ellipse (histogram). When 
several galaxies pass this criterion, the bary center of the group 
is taken. The model (curve) assumes Gaussian positional errors 
and take into account contamination by "chance" counterparts 
6% at center of ellipse) . 

Equation (jij) defines the probability density for a gen- 
uine identification to lie around the IRAS detection. We 
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can compare it to tlie probability of encountering by 
chance an optical object of class c, brighter than mag- 
nitude m between r and r -\- dr: 



Pf (r, TO, c) dr = 2TraxO'YNc{m) r dr, 



(5) 



where Nc{m) is the projected density of optical objects 
with class c brighter than m. This leads to the likelihood 
ratio definition from Wolstencroft et al. ( 1986| ): 



L (r, TO, c) 



Pr{r,m,c) 

Pf{r,m,c) 2TraxcrYNc{m) 



Sutherland & Saunders (1992) stressed that L is inde- 



pendent of the number of optical candidates lying in the 
error ellipse. Some extra information can be gained from 
the presence or absence of these other candidates if one 
assumes that the true optical counterpart to the infrared 
source is unique. However, this hypothesis is dangerous 
for IRAS galaxies, because of the non-negligible fraction 
of m ultiple systems optically identified (e.g. Gallimore & 
Keel 1993| ). Besides, the artificial and real optical galaxy 
samples do not have the same density of "background" 
galaxies, as the first one comes from an infrared flux- 
limited sample: any estimator that combines information 
from the other optical detections is likely to behave differ- 
ently in both cases. We therefore conformed to the "stan- 
dard" likelihood ratio L as an individual identification es- 
timator for optical candidates. 

For every IRAS detection lying within the area of our 
optical catalog, likelihood ratios Li associated to every 
potential optical candidate were obtained in the follow- 
ing way. First, we estimated the Nc(Bj) separately for 
stars (c — S) and galaxies (c — G) from the number 
counts relative to the current plate. Then we defined a 
circular zone with a 3' radius around each IRAS posi- 
tion. This corresponds to a large confidence interval for r: 
at least 6cr, or a diameter 3 times larger than the IRAS 
beam FWHM. Finally, the Li were computed using (||) 
for all optical detections within the circular zone. The 
same procedure was applied to the artificial optical sam- 
ple around detections in the simulated 60 /im images, but 
this time using Ncim) = dexp0.5(Bj — 15.0), which pro- 
vides a good fit to our integrated galaxy number counts 
per square degree (paper I). The sizes of the "artificial" 
error ellipses were computed with formulae (II. F. 5) and 
(II.F.6) in Moshir et al. ( |1992| ), assuming the number of 
scans to be (N) = 20, and a background noise of 30.2 mJy. 

5.4- Identification rate 

What real fraction of infrared sources possess one or more 
optical counterparts? Using the individual likelihood ra- 
tios Li, we can define a global estimator £ for each IRAS 
detection, from which we will derive an identification ra- 
tio. Here again, we need to be able to compare statistics 



done with different densities of contaminating sources. For 
this reasons, we have chosen 



C = sup (Li 



(7) 



though it might certainly be less effective for multiple 
sources than, e.g. C = J2i^i- presence of contam- 

ination by "chance" optical objects, the observed proba- 
bility for C [as defined by (0)] to exceed some threshold 
Cq for an infrared detection reads: 



(6) p{C> Co) = Pr{C > Co) (1 -Pf{C > Co))+Pf{C > Co),{8) 



where Pr{C > Co) and Pf{C > Co) are the probabilities 
to have C > Co because of real and fortuitous optical 
identifications, respectively. Thus 



Pr = 



P-Pf 



1 



Pf 



(9) 



Pf can be measured by putting error ellipses at random 
positions in both real and artificial catalogs, which then 
enables us to compare the real identification rates Pr- 




0.1 1 

/^(60/xm) (Jy) 

Fig. 4. Identification rate (corrected for the effects of spurious 
identifications) of IRAS sources with a likelihood threshold of 
5. The lines are the Monte-Carlo models, as in Fig. 0. 

Figure ^ shows Pr as a function of the 60 fj,m flux, for 
Co — 5 (which proves to be a good compromise between 
completeness and reliability) . Although they do not enable 
us to choose between the models, these statistics seem 
to indicate that the faintest infrared sources cannot be 
attributed to a large number of unexpected "noise peaks" . 
Another interesting result is the rather low identification 
rate predicted by simulations at high S/N: about 80%. 
This in good agreement with the observations, and the « 
90% identification rates report ed for the FSC (Sutherland 
et al. 1991, Wolstencroft et al. 1992), once the differences 



in the sizes of error ellipses and the addition of chance 
identifications have been accounted for. An examination 
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of the simulated FSS images reveals that the unidentified 
sources are produced by confusion with close neighbours in 
projection. Because of the large width of the IRAS beam, 
the barycenter of some detections can be shifted by a few 
pixels, i.e. 5a or more at high S/N, adding a a long tail to 
the Gaussian error distribution. For these detections, the 
measured likelihood ratio of the optical counterpart (s) is 
thus very close to zero. We conclude that this phenomenon 
might be responsible of most of the high S/N unidentified 
point-sources in the FSC. Note that cirrus sources do also 
constitute potential deflection factors. 

5.5. Optical colour distribution 

The colour distribution for all identifications with £ > 10 
(roughly a 10% contamination) is shown in Fig.||, together 
with the average colour tracks produced by the differen- 
tial fc-corrections computed with the Pence ( [1976[ ) SEDs. 
For clarity, only the no-evolution model is plotted. The 
(3 = 3,4 models yields tracks very similar to those of the 
no-evolution model, which is normal as Eq. (^ provides 
very little evolution in optical luminosity for luminous FIR 
sources. The dispersion in colours appears to be substan- 
tially smaller than for an optically limited sample, and as 
can be seen, the bulk of the IRAS population follows quite 
closely the Sbc track, supporting our choice in §5.2.2 to 



adopt an Sbc type for the average optical fc-correction. 



I 1 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




- 
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12 14 16 18 20 22 
B, 



Fig. 5. Colour-magnitude diagram of all optical counterparts 
to IRAS sources with £, > 10, in the photographic B,j and 
R bands (points). This high threshold in likelihood limits the 
proportion of fortuitous sources to about 10%. Solid curves are 
the average colour tracks predicted by the non-evolving model 
described in the text, assuming different fc-corrections. From 
top to bottom: E/SO, Sab, Sbc {heavy line), Scd and Sdm. For 
comparison, the dashed line indicates the average colour of our 
complete optical galaxy sample as a function of magnitude. 



5.6. Magnitude distribution 

We now turn to comparing the magnitude distribution 
with model predictions. Once again, we have to subtract 
the contribution from fortuitous detections. One can de- 
compose the total observed number of identifications per 
magnitude bin at magnitude m, n{m), as 

n{ra) = (1 — pf)nr{m) + (1 — pr)nf{m) + ndm), (10) 

where nr{m) is the number of real counterparts and nf{m) 
the number of fortuitous optical galaxies in the same mag- 
nitude interval. nc{m) represents the contribution from 
IRAS sources for which both real and fortuitous identifi- 
cations can occur (with probability p^P/ ) • We don't know 
Tic (to), but we can safely bracket it between the two ex- 
treme cases: 



< nc{m) < pfnr{m) + prnf{m). 



(11) 



The lower bound corresponds to the case where the mag- 
nitude bin is totally depleted by the presence of the two 
"competitors", and the higher bound corresponds to the 
case where there is no depletion at all. The true value is 
quite probably intermediate. Inserting nc{m) in ( pO| ) we 
get 

nim) - (jE^) rif 
n{m) — nf{m) < nr{m) < ■ — . (12) 

The magnitude distribution of IRAS sources with £ > 5 
is plotted in Fig. H, together with model predictions, for 
two 60 /im flux domains. As expected from the evolution 
models, excess sources are found on the dim side of the 
histogram. A very good agreement is found for 3 < Q < 4, 
but like the identification rate statistics, they do not put 
further constraints to the exact amount of evolution. 

6. Summary and comparison with other works 

Although they cannot be considered as really independent 
one from each other, the results of the four statistical tests 
described above (background fluctuations, number counts, 
identification rate and magnitude distribution) do provide 
a picture remarkably consistent with the phenomenologi- 
cal L* oc (1 -I- z)^ evolution model, with 3 < Q < 4. The 
most stringent test is given by the number counts with 
Q — 3.2 ± 0.2 ± 0.3 down to 0.15 Jy, with some indication 
of a stronger evolution below 0.2 Jy. However, the fluc- 
tuation analysis indicates that Q cannot increase much 
beyond Q = 4 below 0.1 Jy. As indicated earlier, a pure 
density evolution yields a similarly good fit to the data, 
with a comoving density evolving as oc (1 -I- z)^-^ ^^■^^'^■'^ 
down to 0.15 Jy. 

These values are in good agreement with the value 
found on shallower surveys by Saunders et al. ( 1990 ) who 
found Q = 3.1±1.0 — altho ugh according to some au- 
thors, e.g. Fisher et al. 1992 or Oliver et al. 1995 this 
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Fig. 7. Corrected VFSS number counts compared with previ- 
ous studies and simple no-evolution (NE) and L* oc (1 -I- 2)'' 
models at 60 nm. The correction was done by subtract- 
ing to each bin (/i/(60 fj,m)< 0.3 Jy) the bias measured in 
Monte-Carlo simulations, with the Q = 3.5 model. This bias 
depends weakly on the exact value of Q; for 3 < Q < 4, the 
typical uncertainty is ~ 5%, and has been included in the error 
bars. Note that the other counts have not been corrected for the 
Eddington bias, and might be overestimated in their faintest 
bins. The grey area indicates approximately the la confidence 
domain for number counts in the context of a L* ex (1 + z)'' 
pure luminosity evolution model, drawn from the background 
fluctuation analysis (§[l.l|). 



Fig. 6. Magnitude distribution of the optical counterparts 
to IRAS sources in the photographic B,j band in two in- 
frared flux domains. Up: 0.2 < /^(60/im) < 1 Jy, down: 
0.11 < /y(60/im) < 0.2 Jy. The two histograms are the lower 
and higher limits to the observations induced by the correction 
of spurious identifications (see text) , but do not include Poisso- 
nian uncertainties. The lines are the models as in Fig. |l]. Note 
that the optical catalogue is expected to become incomplete 
for Bj > 20 and Bj < 16. 



number might largely be affected by biases — , and Oliver 
et al. ( |1995D with g = 3.3 ± 0.8. However they do a prion 
rule out the extreme evolution discovered by Gregorich 
et al. ( |1995| ) in their IRAS number counts. In particular, 
their high density of faint sources is incompatible with 
our background fluctuation analysis. A check done on the 



Gregorich et al. (1995) fields reveals that they contain cir- 



rus, which might produce a large fraction of false detec- 
tions, and should also contribute significantly to the to- 
tal background noise and its associated Eddington bias. 
Strong cirrus emission also affects the Hacking & Houck 
( |l987|) deep field, but the dispersion in their scanning an- 
gles smoothes considerably edge-effects on cirrus struc- 
tures. This might explain why their counts are in better 
agreement with ours. Our counts and those from previous 
studies are summarized in Fig. |^. 



7. Discussion 



7.1. The nature of the faint source excess 



Our statistics probes mainly the evolution of infrared lu- 
minous galaxies (Lgo > L*). The good fits given to the 
data by the L* oc (1 -|- z)^-^ model, and the large excess 
of faint optical counterparts indicate that luminous IRAS 
galaxies should be « 5-7 times more numerous at z « 0.3 
than at the present time. What are these galaxies? Ac- 
cording to our Monte-Carlo model, the optical counter- 
parts with 20 < Bj < 21 should have an average 60 /im 
luminosity Lgo ~ 5.10^^h~^LQ and should be at a red- 
shift z « 0.3, in agreement with their rather red colours 
Bj — Kp ~ 1-6, if their optical continuum is close to 
that of local Sbc's. These are therefore potential Ultra- 
Luminous InfraRed (ULIRGs) candidates. This picture is 
consistent with the redshift distribution of the op tically 
faint ULIRGs surveyed by Clements et al. ( 1996a ). Ob- 
jects with brighter optical counterparts (Bj < 20) have 
naturally a lower FIR excess, but most are still redder 
than the average population in the same magnitude do- 
main. They are very probably similar to the population 
of starburst spirals with w M* optical luminosities iden- 
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tified in the sub-milliJansky radio survey^ (Bcnn ct al 



1993, Windhorst ct al. 1995), ahhough at smaller redshifts 
(z « 0.1-0.2). Hence it seems unlikely that a large frac- 
tion of the detected excess of faint 60 fim sources at the 
100 mjy level might be caused by dust-enriched primeval 



galaxies, as suggested by Franceschini et al. (1994). 

Wang (|19914 |l991lj ) proposed a pure FIR-luminosity 
evolution scenario based on the evolution of dust content 
in normal galaxies. Although some variants of this model 
might be able to predict a significant increase of FIR lumi- 
nosity in galaxies over the few past Gyr, they rely on the 
assumption that emitters are optically thin at blue/UV 
wavelengths. This is certainly not the case for most of FIR- 
luminous galaxies (the sole population of galaxies that can 
be detected by IRAS at significant redshifts) . As these ob- 
jects already have Lfir > , adding more dust should 
not be very effective in increasing their FIR luminosities 



The z < 1 redshift distributions obtained from spectro- 
scopic surveys are however apparently in good agreement 
with no-evolution models, which suggests a diminution or 
a fading with time of low and moderate luminosity mem- 
bers of the blue galaxy populationci since z w 0.5 (Broad- 
hurst et al. |1984 Lilly |1993| , Ellis et al. |1996| , Driver et al. 



1995 ). Various scenarios to explain this observational phe- 



nomenon have been proposed, including a delayed forma- 
tion epoch o f dwa rf galaxies (Babul & Rees 1992, Babul 
& Ferguson 199(]| ), decimation through intense merging 
(Guid erdoni & Rocca-Volmerange 1991, Broadhurst et al. 
1992), or more prosaically, biases in the determination 



of the local luminosity function (Davies 199C, McGaugh 



(e.g. Belfort et al. 1987). Therefore a moderate change in 
the dust content of galaxies since z « 0.3 should not be a 
major contributor to the evolution observed yet for IRAS 
galaxies. 

The evolution rate of luminous IRAS galaxies is in re- 
markable agreement with those of QSOs (e.g. Boyle et al. 
or X-ray selected Active Galactic Nuclei (AGN) 



1994 Gronwall & Koo |1995D . 

Now, strong evolution among the spiral population, 
apparently seen at FIR and radio wavelengths, is a tempt- 
ing argument for hypotheses advocating recent evolution 
of the optical LF. Evolution in FIR is detected in red- 



shift survey s of IRAS samples (Saunders et al. [1990| , Oliver 
et al. 1995) and could be explained, for example, through 
an increasing rate with z of galaxy-galaxy interactions, 
which are known to be efficient triggers of starburst ac- 
tivity (e.g. Larson & Tinsley 1978| ), or AGNs (Sanders 
et al. 1988] ). There are indeed both theoretical (Toomre 



e.g. Maccacaro et al. 1991), which brings up the link be- |1974 Carlberg£990|) and observat ional (Zepf & Koo |1989 



tween AGN and ULIRGs, as most of the latter may har- 



bour AGNs (e.g. Sanders et al. 198S). An important ques 



tion is therefore: does the evolution concern the global 
FIR-luminous galaxy population, or is it limited to "mon- 
sters" containing active nuclei? The magnitude distribu- 
tion in our faintest 60 ^m flux bin (Fig. ^ seems to in- 
dicate that most of the sources found in excess are not 
ULIRGs, and cons equently, according to local statistics 
(Leech et al. 1989 ), are less likely to contain an AGN. 
A stronger argument against a sample dominated by 
AGNs may come from spectroscopic identifications of sub- 



Burkey et al. 1994, Carlberg et al. 1994) evidences for such 
a rapid increase with redshift of the merging rate, at least 
to z « 0.7. A few attempts have been made to link a 
strong, starburst-driven, evolution of the FIR luminosity 



function to that of the optical (Lonsdale & Harmon 1991 
Lonsdale & Chokshi 1993| , Pearson & Rowan-Robinson 
19961 ). They involved simple proportionality between the 



milliJansky radio-sources (Benn et al. 1993), which, as we 
saw, can be compared to our sample, and has proven to 
contain a minority of AGN-like spectra. Extreme caution 
should however be made in interpreting these spectra, as 
optical signatures of a central engine might be hidden by 
the large amounts of dust present in the nuclear regions 



of such objects (see e.g. Dennefeld 1996). 



7.2. The link with optical evolution 

Deep number counts conducted at visible wavelengths 
have revealed a large excess of faint galaxies over the no- 
evolution predictions from qq > cosmological models 
(Tyson |1984 CoUess c t al. |1990| , Lil ly et al. |1991| , Met- 



calfe et al. 1991, 1995, Small et al. 1995| ), as well as a 



blueing of the global population with fainter magnitudes. 



optical and FIR evolution (either in density or luminos- 
ity) and therefore predicted a large increase in the comov- 
ing space density of optically bright (M < M*) galaxies 
with redshift, which is unfortunately not seen at z < 0.5. 
The same problem arises with stellar population synthesis 
models applied to starburst-driven evolution (e.g. Carl- 
berg & Chariot |l99^ . 

But this apparent contradiction could be alleviated by 
remarking that the possibly large increase of blue con- 
tinuum emission resulting from starburst activity is often 
largely hidden at blue/UV wavelengths in FIR luminous 
galaxies. It is a fact that the complex role of dust, and 
in particular its geometry with respect to the stars, has 
been ignored in many models trying to reproduce number 
counts with starburst-evolution. In massive spirals, star- 
formation induced by merging is generally concentrated in 
the central regions of the galaxy (e.g. Hummel 1981, Con- 
don et al. 1982, Bushouse 1987). Such a feature is likely 
to originate from radial inflow of disk gas triggered by 
the interaction, and is well reproduced in numerical sim- 



ulations (e.g. Hernquist 1989, Minhos & Hernquist 1994) 



Incidentally, the evolution of sub-mjy radio sources seems 
comp atible with that of IRAS galaxies (Rowan- Robinson et al. 



comp a 
199^ ). 



^ For simplicity, we shall refer to these objects as "faint blue 
galaxies" . 



E. Bertin, M. Dennefeld, and M. Moshir: faint IRAS sources 



11 



The associated circumnuclear starburst (which may sur- 
round an AGN) is usuaUy heavily reddened in the optical, 
and does not increase so much the total blue/UV light 
output of the galaxy: ULIRGs are rarely "ultraluminous" 
at visible wavelengths (as reflected by relation |^) . On the 
contrary, starbursts in field dwarf galaxies, dominated by 
gas-ric h, du st-deficient late ty pes (e .g. Van den Bergh & 
Pierce 1991 , Wang & Heckman 1996| ) typically sho w up i n 
optically prominent giant H II regions (e.g. Huchra |l987|) , 
and are often found to dominate the total blue/UV emis- 
sion. Hence, statistically, starbursts observed in both types 
of galaxies possess very different observable signatures. 
Therefore the hypothesis that starburst galaxies which 
dominate IRAS faint number counts are the same as those 
that make the excess seen in the optical at intermediate 
redshifts (Pearson & Rowan- Robinson 1996) is not con- 
vincing. 

There is nevertheless evidence from optical redshift 
surveys that the galaxies with M > M* observed in excess 
at z ^ 0.5 are actively star-forming galaxies (Ellis 1996| , 
and references therein). Recent studies have estimated 
their evolution in terms of global brightening at blue wave- 
lengths between now and z « 0.3 — 5 to ~ 1 — 1.5 mag 



ing properties barely distinguishable, within measurement 
errors, from their common parent population of late-type 



galaxies (e.g. Infante & Pritchet 1995, Mann et al. 1996 



(Driver et al. 1996, Rix et al. 1996), similar to what we 
inferred for IRAS galaxies in FIR. Is the evolution of 
FIR-luminous galaxies closely related to that of faint blue 
galaxies? In the hypothesis of a recent evolution dom- 
inated by starburst processes in both populations, this 
question is essentially linked to the existence of a com- 
mon triggering mechanism, in which case we may consider 
IRAS galaxies as massive, dusty versions of the faint blue 
galaxies. Although interaction/merging seems to be the 
rule among distant IRAS galaxies, at least for the most 
luminous members (Clements et al. 1996b , 1996c ), the sit- 
uation is much less clear for faint blue galaxies. Deep Huh- 
ble Space Telescope images present visual evidence for a 
high proportion of interacting/merging objects, increasing 
steeply with magnitude (Driver et al. 1995, van den Bergh 
et al. 1996). Nevertheless, many of these sources may well 
be luminous galaxies at much higher redshifts, or in some 
cases, knots of star-formation observed in single objects. 
Besides, optical observations of nearby dwarf galaxies un- 
dergoing starbursts reveal that a large fraction are ap- 



parently isolated objects (Campos- Aguilar & Moles 1991 



Telles & Terlevich 1995), which argues against interac- 



tions as a unique triggering mechanism of starburst activ- 
ity at recent epochs. Still, one cannot exclude the hypoth- 
esis of interactions or collisions with HI clouds, supported 
by 21 cm observations of th e environment of nearby HII 
galaxies (Taylor et al. 1995 ). 

Another test of the interaction scenario might be pro- 
vided by studying the spatial (or angular) two-point cor- 
relation function. One could expect the distribution of 
sources triggered by interactions to exhibit some differ- 
ences with respect to the global (late- type) density field. In 
fact both IRAS and blue galaxies happen to have cluster- 



Oliver et al. 1996), although this may not be true at very 



small scales (Carlberg et al. 1994, Infante et al. 199£) 



We shall therefore conclude this discussion by point- 
ing out that a simple and appealing scenario, in which 
the strong, recent evolution among late-type galaxies at 
moderate redshifts is related to an increase of the merg- 
ing/interaction rate with lookback time, appears qualita- 
tively compatible with the observations. In this scenario, 
FIR-luminous galaxies would then represent the massive, 
dusty counterparts of faint starbursting blue galaxies. 
More data are clearly necessary to test this hypothesis on 
the observational side, under others to follow the evolu- 
tion of the intermediate luminosity population in the FIR. 
This should be possible with the ISO {Infrared Space Ob- 
servatory) satellite. It would also allow to assess more ac- 
curately the evidence for strong evolution at low redshifts 
in the FIR, deduced so far only from IRAS measurements. 

8. Conclusions 

We have analysed a deep homogeneous 60 fj,m subsample 
of galaxies, at the lowest flux limit reached by the IRAS 
Faint Source Survey, yielding the following results: 

1. Detection and measurement biases were assessed 
through Monte-Carlo simulations, which prove to repro- 
duce satisfactorily most observed features of detected 
IRAS sources at 60 fim. In particular, the rather high pro- 
portion of unidentified sources in the FSC (> 10%) can 
totally be accounted for by defiections caused by neigh- 
bouring IR sources. 

2. Background fluctuation analysis and number counts 
provide evidence for strong evolution in FIR luminosity 
cx (1 + z)3-2±o.2±o.3 density oc (1 -f z)6-0 ±0-5±o.7^ 

agreement with previous studies. 

3. Sources in excess are genuine and are generally as- 
sociated with faint, relatively red, optical counterparts 
which we interpret as being M* massive starbursting 
galaxies at redshifts ^0.1. 

4. Photometric properties of distant IRAS galaxies 
suggest that the advocated starbursts (or AGNs) are 
considerably hidden by dust at UV/visible wavelengths, 
which would explain why no large number of related op- 
tically luminous objects is showing up in optical surveys 
at z < 0.5. 

5. Faint IRAS galaxies in excess are unlikely to be 
the ones making the excess in number counts at blue 
wavelengths and z < 0.5, but this does not exclude that 
both populations follow the same evolution mechanism 
(i.e. interaction- induced starbursts). 
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